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Like ants, some microorganisms are known to leave trails on surfaces to communicate. We explore 
how trail-mediated self-interaction could affect the behavior of individual microorganisms when 
diffusive spreading of the trail is negligible on the timescale of the microorganism using a simple 
phenomenological model for an actively moving particle and a finite-width trail. The effective 
dynamics of each microorganism takes on the form of a stochastic integral equation with the trail 
interaction appearing in the form of short-term memory. For moderate coupling strength below an 
emergent critical value, the dynamics exhibits effective diffusion in both orientation and position 
after a phase of superdiffusive reorientation. We report experimental verihcation of a seemingly 
counterintuitive perpendicular alignment mechanism that emerges from the model. 

PACS numbers: 87.18.Gh, 87.17.Jj, 87.10.Ca 


For many animals and microorganisms it is advanta¬ 
geous to know where their companions or they themselves 
have been m- To this end, many creatures leave trails 
of some characteristic substance. A well studied example 
is the pheromone trails of ants [SI [10], which allows them 
to collect food efficiently. Single cell organisms are also 
known to leave trails [iii[ia. It is believed that the 
trails help them form aggregates in sparse populations 
[3 nuns], whereas in denser populations, colonies could 
also result from the combined effect of surface-bound 
motility and excluded volume interactions [13 [ig. For 
bacteria, these trails are often subsumed as exopolysac¬ 
charides (EPS) [13 but may also contain proteins [T7|. To 
be evolutionarily favorable, the (energetic) costs incurred 
by trail formation should balance the advantages gained 
through this form of communication. 

Chemotaxis is commonly mediated by rapidly diffusing 
signalling molecules [3 but, more generally, cell-cell sig¬ 
nalling can also be mediated by trails of macromolecules 
that diffuse much more slowly than the microorganism 
or form stable gels [ISKIHI- Ghemically-mediated inter¬ 
actions between bacteria or eukaryotic cells [3 fT ^]]23 as 
well as artificial active colloids [25l - [27] lead to a variety 
of collective phenomena including collapse, pattern for¬ 
mation, alignment and oscillations. Auto-chemotactic 
effects have been studied in the context of swimming bac¬ 
teria [28l [29] and Dictyostelium cells [30] . While much 
is known about the chemotactic machinery in bacteria 
[23 [S3 and eukaryotic cells [22], relatively little is known 
about trail-mediated interactions. 

Whereas ants have antennae that are spatially well sep¬ 
arated from their pheromone glands [33] , such a clear sep¬ 
aration is difficult for single-celled organisms [nilS3l35]- 
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In addition to sensing the trails left by other individuals, 
microorganisms are also immediately affected by their own 
trails. This suggests that trail-mediated self-interaction 
could play a significant role in the behavior of microorgan¬ 
isms. For example, by providing a mechanism to tune the 
effective translational and orientational diffusivities, or 
by creating distinct modes of motility, and consequently, 
the search strategy. 

In this Letter, we discuss a simple but generic model of 
a microorganism experiencing trail-mediated interactions. 
Focusing on a persistent EPS trail with vanishing diffu- 
sivity but taking its finite width explicitly into account, 
we focus on the immediate self-interaction that previously 
had to be excluded a priori by an ad hoc refractory period 
[22l [36] . While previous work has mostly considered a 
concentration dependent speed [JO] EH El] , a coupling to 
the orientation arises naturally pn \T7l 127] . We find that 
the self-trail interaction modifies the translational and 
orientational motion of the microorganisms and renormal¬ 
ize the corresponding diffusion coefficients, at the longest 
time scale (see Fig. [^. 

Microscopic Model. — We consider a single particle of 
width 2R whose state at time t is defined by its position 
r{t) and orientation n(t) = (cos (/j, sin </?). We model the 
dynamics by prescribing a fixed characteristic speed, vq, 
for the particle, namely 

dtr{t) = von{t). (la) 

The motion will typically be generated via the cooperation 
of a number of molecular motors, whether it is realized 
by the retraction of pili [42l [43] , the extrusion of slime 
|44j . or any other mechanism. This implies significant 
noise in the propulsion force, and consequentially, a finite 
directional persistence. For the simplest, trail-free case, 
we model the orientational dynamics as a purely diffusive 
process, dt(p{t) = C(^), where ^{t) is a Gaussian random 
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FIG. 1. (color online). Sample tr ajectories generated by the 
effective dynamics, Eqs. ( la|lc|ld I, over a period of time lO^r 
(color coded) for: no interaction with the trail, fir = 0 (a), 
weak interaction, fir = 1.2 (b), strong interaction, fir = 1.85 

(c) , close to the localization transition, and above it, fir = 2.15 

(d) . The rotational diffusivity is set to Il(!r = 10“^, and 2r 
is the trail crossing time. Panel (e) shows a magnification 
of the end of trail b with the trail field, ip{r,t), of width 2R 
in green and the current orientation of the microorganism 
(ellipse), n. (f) Schematic depiction of a microscopic model 
system such as P. aeruginosa that uses pili for motility and 
sensing, (g) A schematic for the definition of AS, which is the 
angle between the current body orientation and the bacterial 
trajectory (dotted line). 


variable obeying = 2D^6{t — t') and 11° is the 

microscopic rotational diffusion coefficient controlling the 
persistence time This trail-free model displays a 

translational mean-square displacement (MSD) 5r‘^{t) = 
([r(t) — r(0)]^) that crosses over from ballistic, 6r‘^{t) = 
for, D°t «C 1, to diffusive behavior, 6r‘^{t) = 4Zl°t, 
for D^t 1 where Zl° = Vq/{2D°). Fluctuations in vq 
could also be taken into account in a straightforward 
generalization |^ . 

The trail excreted from the microorganism can be char¬ 
acterized by the density profile '0(r, t) that satisfies the dif¬ 
fusion equation dtip{r, t) — 'Dp\7^tp{r, t) = kSj^ (i— r{t)), 
where k is the deposition rate and {r — r(t)) is a “reg¬ 
ularized delta function” that accounts for the finite size 
i?, and traces its position [normalized as / (Pr5\{r) = 1]. 
Setting Up = 0, we find for the trail profile at time t and 


position X as 

= k f dt'5\{x — . (lb) 

Jo 

We choose a rectangular source, 5\{r) = 0(i?^ — r^)/7ri?^, 
where 0(a:) denotes the Heaviside step function and 
r = jrj. The trail width 2R defines a microscopic time 
scale r = R/vq, which gives the trail crossing time (see 
Fig. [^). This specific regularization scheme is a good 
representation of the regime in which the characteristic 
diffusion length of the polymeric trail is much smaller 
than the width of the trail, R [46] . 

A generic interaction with the trail couples to gradients 
of the trail field perpendicular to the current orientation 
mi HZ], effectively steering the microorganism toward 
trails by favouring a orientation fi perpendicular to the 
trail, i.e.. 


dMt) = xd±'<Pir{t),t) + ^(t), 


(Ic) 


where i9j_'0 = n±{t) ■ with fly = 

(—siny), cos being the angular unit vector in polar 
coordinates. The sensitivity to the trail is controlled by a 
parameter y. We have provided a microscopic derivation 
of this coupling |37j for a model system of a pili-driven 
bacterium on a substrate (see Fig. S) by assuming a 
generic dependence of the pili surface attachment force 
on the EPS concentration. However, Eq. ( [T^ will be ex¬ 
pected in the continuum limit for any microscopic model 
based on symmetry considerations p7] . 

Effective Dynamics .— We assume that the particle 
trajectory does not bend back on itself immediately (no¬ 
small-loops assumption), and that self-intersections on 
longer times are rare enough to be negligible. 

By making a short time expansion of Eqs. (la) and ([lc| 


to be inserted into Eq. (lb I, one finds a closed equation 


for the head of the trail field = d±ip{r{t),t) [57] . 

The result is a stochastic integral equation 

d±'f{t) = — [ du{T-u)[d^il){t-u)+^{t-u)/x], (Id) 
Jo 

The effective turning rate H = kxT/irR^ increases for 
more intense trails (larger k) and for more sensitive organ¬ 
isms (larger y). The delay t reflects the memory imparted 
by the trail. The closed set of equations ( |la|lc|ld| ) con¬ 
stitute our effective dynamical description of the system. 

Eor the average gradient one finds {d±'f) ~ exp(at) 
where the rate a is given implicitly as the solution of 
X{a) = 0 where X{a) = 1 — ^ [l + ^ (6““"” — 1)]. For 
Ht < 2, a < 0 and Eq. (Id) defines a random process 


with zero mean that leads to a stationary dynamics which 
is time-translation invariant. For Hr > 2 one finds a > 0, 
i.e., the gradient (angular velocity) diverges exponentially 
in time. Implying that the trajectory converges in a 
logarithmic spiral to a localized point. This means there 
is a maximum value of the product of trail deposition rate 
and sensitivity, fcy, that allows steady-state motion. For 
sample trajectories, see Figs. and 
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FIG. 2. (color online). Angular MSD (a), and translational MSD 5r^{t) normalized by the trail width, R, (b) as a 

function of time t for several values of the effective turning rate fir = 0,0.4,0.7,1.0,1.3,1.6,1.9. The microscopic diffusivity 
is set to D^t = 0.1. The inset of (a) shows ‘5(p^{t) for fir — 1.99 demonstrating the intermediate superballistic regime, (c) 
Color coded translational diffusivity D normalized by the trail width R and the trail crossing time r as a function of the control 
parameters. Note the logarithmic axes. Inset: The translational diffusivity for the same range of parameters normalized by the 
trail free (fir = 0) value on a linear scale. The white contour lines are 0.1 apart. 


In the stationary regime, Eq. can be solved in the 
frequency domain, x9±ipiuj) = — l]^(a;). The 

trail-mediated self-interaction thus linearly transforms 
the intrinsic white noise ^ to an effective colored random 
angular velocity, S(a;) = such that dtip{t) = 

E{t). 

Angular & Translational MSD. — The most easily 
accessible quantity in experiments is the translational 
MSD 6r'^{t), which is related to the angular MSD, 
= (bW - via gZ] 

Sr^{t) = 2vl f dt' it- (2) 
Jo 

The angular MSD is a sum of three terms, given in the 
Laplace domain as s^(5(^^(s)/(2Zl°) = 1-|-A(s)-|-A(s). The 
corrections to simple diffusion are given by the two corre¬ 
lation functions A{t) := X {9±'J’{t)^{0)) /(2Z1°) (such that 
A(s) = A“^(s) — 1) and A{t) = x^ {d±tl}{t)d±'tj}{Q)) /D^ = 
2 f* dt'A{f)A{t' + t). The definition of A(s) shows that 
fir is the only relevant control parameter for the orien¬ 
tational MSD and the behavior of is fully 

determined by the analytic structure of A“^(s). 

In the stationary regime, (cf. Fig. gB]) starts 

off diffusively, S(p^{t) = 2D^t for t r and becomes 
asymptotically diffusive again, Sip'^{t) = 2Drt for t ^ 
t/(1 — Dt/2) determined by the smallest pole of A“^(ia;). 
The effective orientational diffusivity 


fl > D^, the intrinsic noise combines a diffusive (oc 
excursion with the ballistic (oc t) reorientation due to the 
self-interaction, leading to Sip{t) oc until later times 
where the stochastic character of the self-interaction be¬ 
comes important and turns the behavior back to diffusion. 

The translational MSD 5r'^(t) always starts ballistically, 
6r‘^{t) = for t <C r and crosses over to diffusive 
behavior 5r'^{t) = ADt for long times t —)■ oo (cf. Fig. [!>)• 
The crossover time, t* will be determined implicitly by 
6(p'^{t*) ^ 1. For the location of the crossover and the 
dependence of the translational diffusivity on the control 
parameters, we have to consider a number of different 
regimes. 

Short Persistence Regime. — When < 1 + 

Dt 12, the crossover happens around fit ^ ^/T+~2AlfD^— 
1 and the asymptotic diffusivity 

D/D° = yfTTD'^j2fle’^°/'^^evic{^jD^j2n), (4) 

is a function of the ratio D^/fl alone. 

Long Persistence Regime. — For sufficiently straight 
trails such that by the time 6ip‘^{t) ~ 1 it is already 
deep in the long time diffusive regime, i.e., D^r <C 2(1 — 
Slr/2)^/[2 — S1 t + (flr)^/2], the crossover happens around 
t ^ l/D^ and the asymptotic diffusivity is given as 


D/D° 


^\l- ^T{l + nT/2)^ ' 

Dr 6 (1 — r2r/2)® 


(5) 


Dr/D'^ = 1 + X 


fir 1 + f2r/2 


(1-L!t/2)2’ 


(3) 


diverges for Dr -x 2 confirming our expectation that the 
trail-mediated self-interaction reduces orientational persis¬ 
tence. The two diffusive regimes are joined by an interme¬ 
diate, superdiffusive regime. Note that the crossover time 
to the asymptotic diffusive regime diverges as Hr —>• 2. 
Close to the limiting value Hr = 2, the crossover is 
given by the super-ballistic law 5gP{t) = 6D^T{t/T)^ for 
r<Ct<Cr/(l — Hr/2). For a fast effective turning rate 


Critical Regime. — Close to the upper limit Hr -x 2 
and for D/r < 1, the crossover happens around tjr ^ 
Xj ^D'^t and the asymptotic diffusivity 

D/D° = r(3/4)(D/r)2/3 bItIZ, (6) 

is a function of D/r alone. Note that the dependence on 
the intrinsic noise, D oc 1/ \/D'^, is significantly weakened 
compared to the trail free case, D'^ oc 1/D/, and that D 
does not vanish as Hr —)• 2. 
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FIG. 3. (color online). Phase diagram of the dynamics of 
the microorganism with trail-mediated self-interaction, as a 
function of the dimensionless turning frequency fir. Inset: 
Zoomed view of Fig. Id. 


Intermediate Regime .— In the rest of the parameter 
space of D^t and fir, no explicit expressions can be given 
and the asymptotic diffusivity D will be a function of 
both control parameters. Numerical result for the effective 
translational diffusivity is presented in Fig. W 

Discussion .— The interaction of a microorganism 
with its own trail effectively introduces a new timescale 
1/n. The trail-mediated self-interaction modulates the 
intrinsic noise in a linear but nontrivial way. While the 
asymptotic dynamics remains diffusive below the critical 
value fir = 2, both for the translational and the orien¬ 
tational degrees of freedom, the (orientational) diffusive 
regime may only be reached on timescales that may be 
much longer than r. This holds in the limit of strong 
interactions where the crossover timescale r/(l — D,t/2) 
becomes increasingly large but also for large effective per¬ 
sistence times 11 Dr. Moreover, the crossover times are 
distinct from the microscopic persistence time l/D^. 

The asymptotic angular diffusivity. Dr, is a function 
of rir that diverges as fir —>■ 2. It is always larger than 
its microscopic value D^, i.e., orientational persistence is 
reduced by the trail. The translational diffusivity, D, on 
the other hand, is always reduced and, in general, depends 
on the parameters Hr and D^t. 

For very strong trail-mediated self-interactions, fir > 2, 
the initial dynamics quickly confines the particle in a re¬ 
gion of size ^ 2R. Here, our assumptions break down and 
the ensuing dynamics will depend on more microscopic de¬ 
tails not resolved in the present model. Nevertheless, we 
can formally identify the behavior of the microorganism 
by a diverging rotational diffusivity. The translational dif¬ 
fusivity incurs a finite jump at the transition point as it is 
determined by the regular short-time part of S(p'^ (t) below 
the transition and by unresolved microscopic details above 
it. We note that the no-small-loops condition, which can 
be written as Sip'^{T) ^ (D^t)'^ < tt^, does not ob¬ 

scure this phase transition, so long as D^t < 1. The 
analysis of a more microscopic model is beyond the scope 
of this contribution. The behavior of the system can 
be summarized in a phase diagram as shown in Fig. 
Our findings are subtly different from the sub-diffusion 





FIG. 4. (color online), (a) The experimentally observed nar¬ 
rowing of the A9 distribution with increasing trail deposition 
(see Fig. for the definition of A0). The experiments were 
carried out with the P. aeruginosa mutants Apsl, which does 
not secrete Psl (light green) and APpsi/PBAD-psl, which se¬ 
cretes Psl in response to the arabinose in the environment. 
For AP psl/P bad-PsI, the arabinose concentration was varied 
between 0% (light blue, low Psl deposition) and 1% (light red, 
high Psl deposition), (b) The corresponding theoretical AO dis¬ 
tributions resulting from the influence of the alignment term. 
The trail coupling strength is varied between Qt = 0.05 (green) 
fir = 0.12 (blue) and fir = 0.35 (red). The distribution has 
been sampled using the time interval of At = 0.2 r, and the 
rotational diffusion has been set to DrT = 0.015. These values 
roughly correspond to the experimental parameters at which 
the distribution has been measured, (c) A comparison be¬ 
tween experimental (dots) and theoretical (lines) distributions 
of AO for the AP psi/Pbad-PsI mutant under 0% (blue) and 
1% arabinose (red). 


observed in Ref. m as a result of temporary trapping of 
bacteria in loops caused by quenched disorder, which is a 
stationary regime. 

The perpendicular alignment strategy might appear 
to be counter-intuitive, but it is supported by recent ex¬ 
perimental evidence. We have investigated this question 
by probing the angle distributions of single bacteria in 
experimentally recorded motion of P. aeruginosa with 
different rates of Psl exopolysaccharide deposition (in a 
similar experiment to Ref. [H]); for methods, see the 
Supplemental Material m- For the angle A9 between the 
trajectory and the body orientation (defined in Fig. [^), 
the experiments show that the distribution narrows down 
with increasing increasing the secretion of Psl trails in via 
a mutant with an arabinose inducible promoter, which 
increases the effective trail deposition strength; see Figs. 
1^ and c. Our numerical simulations of the model predict 
a similar trend, as shown in Fig. & and c, highlighting 
the interplay between the noise and the tendency of per¬ 
pendicular alignment. The trail can have a stabilizing 
effect by reorienting the microorganism towards the inner 
regions in case it reaches the trail boundaries. The result 
is a significantly narrower A6 distribution with increasing 
0, indicating a higher correlation between trajectory and 
microorganism orientation. Further evidence in support 
of the perpendicular alignment scenario can be extracted 
from the collective behavior of the bacteria |52j . 

Our results could have significant biological (as well as 
biophysical m) implications. Regulating the strength of 
the trail-mediated self-interaction may allow microorgan¬ 
isms to decide whether to confine themselves to smaller 
areas and search them more thoroughly or explore larger 
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areas. Interestingly, the effective translational diffusion 
coefficient scales as 1/ \/l^ in the presence of strong 
trail-mediated self-interaction, as compared to 1/D^ in 
the trail-free case. This suggests that trails make the 
microorganism less sensitive to intrinsic variations in the 
orientational noise. 

To conclude, we have shown that a very simple model 
of particle-trail interaction leads to a wealth of nontrivial 
phenomena, including an transition from stationary to 
non-stationary behavior with a diverging orientational 
diffusivity. Our results could shed light on the behavior 
of trail-forming microorganisms, and in particular how 
they can use this “expensive” output to regulate their own 


activity while, simultaneously, providing a communication 
channel with other individuals. Moreover, they could also 
find use in the field of robotics by providing a blue-print 
for designing micro-robots that can tune their search 
strategy via local interactions with their own trails. 
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where f = 

We can then use this to approximate the average total 
force (F) = h) 


FIG. 5. Schematics of pili distribution around a bacterium. 

(a) Asymmetric distribution around the body {N. gonorrhoae). 

(b) Pili concentrated at one pole (P. aeruginosa), (c) Idealized 
model with a symmetric distribution of pili at the poles. 


Vol. 73 (Clarendon Press, Oxford, 1988). 

[48] We used the Fixed Talbot method [32] for the numerical 
Laplace inversion and employed SciPy’s [50] quad method 
for the integral in A{t). 

[49] J. Abate and P. P. Valko, Int. J. Numer. Meth. Engng. 
60, 979 (2004). 

[50] E. Jones, T. Oliphant, P. Peterson, et al, “SciPy: Open 
source scientihc tools for Python,” (2001-), http://www 
scipy.org 

[51] O. Chepizhko and E. Peruani, Phys. Rev. Lett. Ill, 
160604 (2013). 

[52] A. Celimson, K. Zhao, C. Lee, W. T. Kranz, G. C. L. 
Wong, and R. Golestanian, to be published. 


Appendix A: Microscopic derivation of the 
equations of motion 

We regard a bacterium of length i which secretes an 
EPS film and pulls itself forward using typically N pili, 
Fig.0 We will assume that the pili attach to the surface 
in an EPS-dependent way. Better attachment results in 
more effective pulling and therefore we can assume that 
the pulling force of an individual pilus attached at site 
ri is EPS-dependent, / = /(■!/'(ri)). The distribution 
of pili on the surface of bacteria varies among species. 
They may be scattered all over the surface (Eig. la), or 
concentrated in one (Eig. lb), or both poles (Fig. Ic) 
[55] . The general argument given below holds for all 
these pilus configurations but to be concrete, we will 
assume a symmetric distribution of an equal number of 
pili, concentrated at the two poles. The attachment sites 
are assumed to be distributed randomly according to some 
distribution P{ri — Tq; n) = P{rp; h) = P{rp,'d), where 
To is the head (or tail) of a bacterium, is the relative 
coordinate of the pilus tip, and d is the angle between 
n and rp. We assume that the force f^ = fWr,))ep 
will pull along the direction of rp. It is plausible to 
assume that the probability for directions in which pili 
can face is symmetrically distributed around the body 
orientation h. This implies that / epP{rp;n)cPrp || n. 
If this was not the case, the pili would preferentially 
explore the space right/left of the body and there would 
be, to a first approximation, a constant torque turning 


(F) J d'^rpP{rp;h)f{'il;{ro + rp))ep 

«F(V'(ro))n f cos i!)P{'0)d-d 


-LF'(V’(ro))VV'(r'o) 


d rpP{vp^ n)TpGpGp 


The pili force component parallel to the bacterial body 
will leave the orientation of the bacterium unchanged and 
propel the centre of mass of the bacterium. The average 
propelling force will be 

(E||) = h{ft ■ (F)) « nF(V’(T-o))(cosi?) (Al) 

In an overdamped system, the velocity of the microor¬ 
ganism will be proportional to the pulling force dtr = 
/ry (F]|) = voh where /ry is the translational mobility. This 
gives equation (la) with 


vo{ip) = {cos'd) (A2) 


The velocity of the microorganism is in general ip- 
dependent but in case of a large constant contribution of 
F, the ^/>-dependent contribution to the velocity will be 
sub dominant. 

The perpendicular component of the pulling force F± = 
F — n(n • F), on the other hand, will generate a torque 
on the body of the microorganism and here, the EPS 
dependence of the force needs to be taken into account 
even at the lowest order. The average torque is given by 
(t) = {ij2)'h X (F). Using Eq. (Al) we get 


(x) « -F'(V’(ro))[n x \/p;{ro)]{rpSiiP-d) 


(A3) 


and a completely equivalent equation for the torque com¬ 
ing from pili pulling at the other tip of the microorganism. 

On a surface the motion of the microorganism will be 
overdamped and therefore the angular velocity uj will be 
linear to the sum of the torques from both tips 


d'h ~ / X 

— = -nXWPS -^_Ln X (xiiead + Ttail) 

= —x^ X (n X Vip) 


(A4) 


where fi± is the rotational mobility and 

X(V') = tJ-±^F'{'tjj{r)){rp sin^ ’d) (A5) 


which becomes independent of ip if F'{ip) « const. 
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In addition to this deterministic reorientation, we 
should also account for noise due to thermodynamic fluc¬ 
tuations and fluctuations in pili attachment, which result 
in torque fluctuations. For the case that the pili pulling 
direction is correlated only on a very short time scale (i.e. 
the pili attach according to P{rp) much quicker than the 
microorganism moves), we obtain the stochastic equation 


Appendix B: Derivation of Equation (Id) 

From the definition of the trail field we have 
2k /■* 

j dt'[x-r{t')]5{R^-\x-r{t')\^) 
and with a change of variables t' ^ t — t', 


2k r* 

dTt ^ 

— =-Xny.{fixV^)+hxir (A6) X(5(i?2-|r(t)-r(t-t')n (Bl) 

From the twice iterated integral 

pt — T pt — T pU 

This equation can be rewritten in terms of the orientation t) —dun{t)+vo Jq 

angle tp relative to some axis on the surface by defining 

n = (cos(p,sin(/3)^, which gives Eq. (Ic). one finds [r{t) — r{t — r)]^ = UqT^ -I- 0 { t ^) and 


pT pU 

[r{t) — r{t — t)] ■ n±(t) = —vq du dw{xez ■ [n{t — w) x — w),t — w)] + — w)}n±{t — w) ■ h±{t) 

Jo Jo 


An identical iteration shows n±{t — w) = l + 0{w) 

and thus 


[r(t) - r{t - r)] • n^{t) 

= —uo du dw[xdj_ip{t — w) + — w)]+0 {t^) 

Jo Jo 


Using the above two approximations in Eq. (Bl) and 
performing one of the integrals yields Eq. (Id). 



FIG. 6. Angular MSD Sp^it) as a function of time t for 
Ot = 1.94, D^t = 0.1 (black crosses) and a power law fit (red 
line) oc for D°t € [0.05,1]. The inset shows the same 
plot over a larger time window, showing that the fit does not 
Appendix C: Implications on experimental capture the asymptotic behavior, 

measurements of diffusivities 


The MSD curve in Fig. highlights subtleties that 
must be considered when interpreting such measurements 
in terms of a model, and for extracting model parameters. 
The long time diffusivity, is not the bare diffusiv- 
ity as determined by the cellular motility apparatus. 
Moreover, the asymptotic diffusive regime is reached on 
very long time scales only: timescales on the order of 
hours |46j . and that may well be beyond experimental 
reach, and even beyond the life time of microorganisms. 
In that case, it might be tempting to conclude that the 
MSDs show anomalous diffusion asymptotically; as can 
be seen from Fig.[^ an anomalous exponent /? « 2.64 can 
fit the data extremely well. 


Appendix D: Experimental Methods and Sample 
Preparation 

P. aeruginosa PAOl |39j strains Apsl - a strain that 
does not produce Psl - as well as AF psi/Pbad-PsI [101 
- an engineered Psl-inducible strain - were used in this 
study. For the detailed experimental information such as 
culture conditions, flow cell assembly and image capture 
system, please refer to methods in m- An overnight 
bacteria culture in FAB medium |41j supplemented with 
30mM glutamate, was diluted and injected into a flow cell. 
FAB medium with 0.6 mM glutamate was continuously 
pumped through the flow chamber using a syringe pump 














with a flow rate of 3ml/hour at 30° C. Different amounts 
of arabinose were added into the medium to control the 
production of Psl. Bright-field images were taken every 3 
seconds by an EMCCD camera on an Olympus 1X81 mi¬ 


croscope equipped with Zero Drift autofocus system. The 
image size is 67 x 67/im^ (1024 x 1024 pixel^). A typical 
data set has about 14000—20000 frames and contains up 
to 1,000,000 bacteria images. 



